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ABSTRACT 

We investigate the generation of intrinsically asymmetric or one-sided outflows or 
jets from disk accretion onto rotating stars with complex magnetic fields using ax- 
isymmetric (2.5D) magnetohydrodynamic simulations. The intrinsic magnetic field of 
the star is assumed to consist of a superposition of an aligned dipole and an aligned 
quadrupole in different proportions. The star is assumed to be rapidly rotating in the 
sense that the star's magnetosphere is in the propeller regime where strong outflows 
occur. Our simulations show that for conditions where there is a significant quadrupole 
component in addition to the dipole component, then a dominantly one-sided conical 
wind tends to form on the side of the equatorial plane with the larger value of the 
intrinsic axial magnetic field at a given distance. For cases where the quadrupole com- 
ponent is absent or very small, we find that dominantly one-sided outflows also form, 
but the direction of the flow "flip-flops" between upward and downward on a time- 
scale of ^ 30 days for a protostar. The average outflow will thus be symmetrical. In 
the case of a pure quadrupole field we find symmetric outflows in the upward and 
downward directions. 
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1 INTRODUCTION 

There is clear evidence, mainly from Hubble Space Telescope 
(HST) observations, of the asymmetry between the approach- 
ing and receding jets from a number of young stars. The ob- 
jects include the jets in HH 30 (Bacciotti et al. 1999), RW Aur 
(Woitas et al. 2002), TH 28 (Coffey et al. 2004), and LkHc^ 
233 (Pererin & Graham 2007). Specifically the radial speed 
of the approaching jet may differ by a factor of two from that 
of the receding jet. For example, for RW Aur the radial red- 
shifted speed is ~ 100 km/s whereas the blueshifted radial 
speed is 175 km/s. The mass and momentum fluxes are 
also significantly different for the approaching and receding 
jets in a number of cases. Of course, it is possible that the ob- 
served asymmetry of the jets could be due to say differences 
in the gas densities on the two sides of the source. Here, we 
investigate the case of intrinsic asymmetry where the asym- 
metry of outflows is connected with asymmetry of the star's 
magnetic field. 
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There is substantial observational evidence that young 
stars often have complex magnetic fields consisting of dipole, 
quadrupole, and higher order poles possibly misaligned with 
respect to each other and the rotation axis (Donati et al. 
2007a, b; 2008; Jardine et al. 2002). Analysis of mater flow 
around stars with realistic fields have shown that a fraction of 
the star's magnetic field lines are open and may carry outflows 
(e.g., Gregory et al. 2006). 

A number of global 3D MHD simulations have been 
done of disk accretion onto rotating stars with complex mag- 
netic fields. The star's intrinsic field may be a superposi- 
tion of aligned or misaligned dipole and quadrupole fields 
(Long, Romanova, & Lovelace 2007, 2008), or a superposition 
of dipole and octupole field components (Long, Romanova, 
Lamb 2009; Romanova et al. 2009a; Long et al. 2010). These 
simulations were focused on accretion processes. To study the 
outflows from these systems requires a much lower coronal 
density than assumed in these works. Intermittent outflows 
from the disk-magnetosphere boundary have been found in 
axisymmetric simulations in cases where the star has a dipole 
magnetic field and where symmetry about the equatorial 
plane was assumed (e.g., Goodson, Winglee, & Bohm 1997; 
Goodson, Bohm, & Winglee 1999). 
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Section 2 of the paper describes the different aspects of 




the simulation model. Section 3 discusses the results of the 
simulations We will quantify the asymmetry of the velocities 
of the top/bottom jets as well as the asymmetry of the mass 
and momentum fluxes for different quadrupole/dipole field 
strengths, different disk accretion rates, and different stellar 
rotation rates. Section 4 gives the conclusions of this work. 



2 SIMULATION MODEL 

We investigate the formation of one-sided or asymmetric jets 
using a set of axisymmetric numerical simulations. The ar- 
rangement of the problem is similar to that described by 
Ustyugova et al. (2006, hereafter - U06). However, in contrast 
with U06, the simulations were performed in the entire sim- 
ulation region without assuming symmetry about the equa- 
torial plane. Below we describe the main aspects of the sim- 
ulation model and point out the new aspects of the present 
model. 



Figure 1. The magnetic field lines ^(r, z) = const and constant mag- 
netic pressure lines for the case of an aligned dipole and quadrupole 
field where the flux function is ^ = /i^^^/^^ + (3/4)/ig2;r^/i?^, 
where B? = r"^ ^ and /i^z is the dipole moment and /ig is the 
quadrupole moment. Roughly, jiq/ jid is the distance at which dipole 
and quadrupole fields are equal. The funnel flow (ff) and the wind in 
this figure are suggested. The dashed lines are constant values of B^. 

In long-time axisymmetric (2.5D) simulations, long- 
lasting outflows were obtained first in the propeller regime 
where the star spins rapidly (Romanova et al. 2005; Ustyu- 
gova et al. 2006), and subsequently in the general case by Ro- 
manova et al. (2009b; hereafter R09). In the propeller regime 
the outer part of the star's magnetosphere - the magnetopause 
- rotates more rapidly than the Keplerian rate of the accre- 
tion disk (Lovelace, Romanova, & Bisnovatyi-Kogan 1999). 
Simulations show that most of the matter outflow goes into 
conical-shape winds from the inner part of the disk. This wind 
resembles the X-wind model (Shu et al. 1994), but there are 
a number of important differences discussed in R09. 

The present work investigates the nature of outflows 
from a star with a complex magnetic field. As a first step we 
consider superposition of an dipole and a quadrupole field 
components both aligned with the rotation axis of the star 
and of the disk. For such a configuration the magnetic field is 
not symmetric about the equatorial plane. This can give rise to 
one-sided or asymmetric magnetically driven outflows (Wang, 
Sulkanen, & Lovelace 1992). Figure 1 shows the nature of a 
combined vacuum dipole plus quadrupole field components. 
In this case we have sketched an extreme possibility where 
there is a conical wind from the top side of the disk but no 
outflow from the bottom side, only a funnel flow marked ff. 

In this paper we present results out systematic MHD sim- 
ulations of the formation of conical winds for the axisymmet- 
ric dipole/quadrupole field combinations following the ap- 
proach of R09 where conical winds are found to form at the 
disk-magnetosphere boundary. The new aspect of the present 
work is that the outflows are not required to be symmetrical 
about the equatorial plane. 



2 . 1 Basic Equations 

Outside of the disk the flow is described by the equations of 
ideal MHD. Inside the disk the flow is described by the equa- 
tions of viscous, resistive MHD. In an inertial reference frame 
the equations are 

£e + v(pv) = o, (1) 

— - Vx(vxB) + Vx (j7tV X B) = , (3) 

^+V-(p5v) = Q. (4) 

Here, p is the density and ^S* is the specific entropy; v is the 
flow velocity; B is the magnetic field; T is the momentum 
flux-density tensor; Q is the rate of change of entropy per unit 
volume; and g = —{GM/r'^)r is the gravitational accelera- 
tion due to the star, which has mass M. The total mass of the 
disk is assumed negligible compared to M. The plasma is con- 
sidered to be an ideal gas with adiabatic index 7 = 5/3, and 
S = \n{p/p^). We use spherical coordinates [r,6,(j)) with 6 
measured from the symmetry axis. The condition for axisym- 
metry is d/dcj) = 0. The equations in spherical coordinates are 
given in U06. 

The stress tensor T and the treatment of viscosity and 
diffusivity are described in the Appendix of R09. Briefly, both 
the viscosity and the magnetic diffusivity of the disk plasma 
are considered to be due to turbulent fluctuations of the ve- 
locity and the magnetic field. We adopt the standard hypoth- 
esis where the microscopic transport coefficients are replaced 
by turbulent coefficients. We use the a-model of Shakura and 
Sunyaev (1973) where the coefficient of the turbulent kine- 
matic viscosity vt = ayC^/^K, where Cs is the isothermal 
sound speed and QK{r) is the Keplerian angular velocity. 
Similarly, the coefficient of the turbulent magnetic diffusiv- 
ity r]t = aacl/QK' Here, ay and aa are dimensionless coeffi- 
cients which are treated as parameters of the model. 
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Figure 2. Development of an asymmetric jet in the main case ip^d = 10, p.q = 20). The color background shows the matter flux-density. The lines 
are the poloidal field lines. The simulations are shown at different times t, which is measured in periods of rotation of the disk at r = 1. 




Figure 3. Enlarged snapshot of the outflow in the main case ip^d = 10, pq = 20) att = 50. The color background shows the matter flux-density 
and the lines are the poloidal field lines. The vectors show the poloidal velocity. 
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Figure 4. The color background shows rB^ (r, z) att = 50, which is proportional to the poloidal current flow through a circular disk of radius r 
at a distance z. The lines are poloidal field lines and the arrows show the poloidal velocity. For this case /ic^ = 10 and fiq = 20. 




2.2 Magnetic field of star 

We consider the superposition of aligned dipole and 
quadrupole field components, 

B = Bd + B, , (5) 
where 

_ 3/id(/Xd • r)r _ i^dfid 
' - |r|3 |r|3 ' 

^ 3/ig(5(/Xg • - l)r _ 3/ig(Ag • r)Ag 
4|r|4 2|r|4 

(6) 

Here, Bd and Bg are the magnetic fields of the dipole and 
quadrupole components with /x^ and /ig the magnetic mo- 
ments. Also, r, fid and fiq are unit vectors for the posi- 
tion and the direction of the dipole and quadrupole mo- 
ments, respectively. For the considered conditions /Id and 
jlq are in the z— direction. The combined dipole/quadrupole 
field can be expressed in terms of the flux function ^ = 
fidr'^/R^ + {3/4:)fiqzr'^ /R^, where the field lines correspond 
to ^(r, = const. Here, we briefly use cylindrical coordi- 
nates (r, 0, z), with = ^ -{l/r)d'^/dz, and 
B, = {l/r)d^/dr. 

2.3 Reference Units 

The MHD equations are solved in dimensionless form so that 
the results can be readily applied to different accreting stars 



(see §7). We take the reference mass Mo to be the mass 
M of the star. The reference radius is taken to be twice 
the radius of the star, Ro — 2 x R^. The reference veloc- 
ity is vo = [GM / RqY^'^ , and the reference angular velocity 
— 1/to. We measure time in units of Po = 27rto, which 
is the Keplerian rotation period of the disk at r = i?o . In the 
plots we use the dimensionless time T — t/Po. 

The dimensionless dipole and quadrupole magnetic mo- 
ments are 

~ = ~ = frj^. 

where Bo is the reference magnetic field. Taking into ac- 
count that fid = BdirRl = BoRojld, we find Bq = 
Bdi.{Rir/Ro)^ /fid, where is the equatorial dipole mag- 
netic field strength on the surface of the star. 

The reference density is taken to be po = Bq/vq. The 
reference pressure is po = Bq. The reference temperature is 
To = Po /^po = vq /IZ, where 71 is the gas constant. The ref- 
erence accretion rate is Mq = pqvqRq, The reference energy 
flux is ^0 = Mo^^o- The reference angular momentum flux is 
Lo — MoVqRq. 

The reference units are defined in such a way that the 
dimensionless MHD equations have the same form as the di- 
mensional ones, equations (l)-(4) (for such dimensionaUza- 
tion we put GM — 1 and IZ — I). Table 1 shows exam- 
ples of reference variables for different stars. We solve the 
MHD equations (l)-(4) using normalized variables: p — p/ po, 
V ^ v/vo, Bd ^ Bd/Bo, Bq = Bq/Bo etc. Most of the plots 
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Figure 6. The figure shows the initial poloidal field lines for a fixed dipole moment /x^ = 10 and different quadrupolar moments, fiq 
0, 20, 30, 40. The field lines have the same value of the flux function in all plots. 




Figure 7. Outflows for different quadrupole strengths /^q but the same dipole component /j.^ 
matter flux-density distribution. The lines are poloidal magnetic field lines. 



10 at t = 50. The color background shows the 



show the normalized variables (with the tildes implicit). To 
obtain dimensional values one needs to multiply values from 
the plots by the corresponding reference values from Table 1. 



2.4 Initial and Boundary Conditions 

The initial and boundary conditions are the same as those 
used in U06. Here, we summarize these conditions. 

Initial Conditions. A star of mass M is located at the 
origin of the coordinate system. A cold disk and hot corona 
are initialized in the simulation region. The disk is of low- 
temperature Td and high-density pd. The corona is of high- 
temperature Tc^Td, and low-density pc <^ pd and it fills all 
other space but the disk. The disk extends inward to a radius 
= 5 and rotates with Keplerian angular velocity oj ^ Qk- 
In reality, it is slightly sub-Keplerian, Q{0 = 7r/2) = kQr 
(/^ = 1 — 0.003), due to which the density and pressure de- 
crease towards the periphery. Initially, at any cylindrical ra- 
dius r from the rotation axis, we rotate the corona and the 
disk at the same angular rate. This avoids a jump discon- 
tinuity of the angular velocity of the plasma at the bound- 
ary between the disk and the corona. Inside the cylinder 
r < rd, the matter of the corona rotates rigidly with angu- 
lar velocity Q{rd) = K{GM/r^y^'^. For a gradual start-up we 
change the angular velocity of the star from its initial value 
Q{rd) = 5"^/^ ^ 0.09 ird = 5) to a final value of 11* = 1 over 
the course of three Keplerian rotation periods at r = 1. 

In most simulation runs we fix the dipole moment of 
the star as pd = 10 and vary the quadrupole moment as 



pq = 0, 1, 10, 20, 30, 40, 60. We also have test cases of pd = 
and pq = 20, 60. The angular velocity of the star in the pro- 
peller regime is = 1 and this corresponds to the corotation 
radius of rcor = 1. The initial density in the disk at the fidu- 
cial point (at the inner edge) is pd — 1, initial density in the 
corona pc — 0.0003. The gas in the corona is hot with initial 
temperature Tc — I and the disk is cold with temperature 
Td = {pc/ Pd)Tc = 3 X 10~^. There is initial pressure equilib- 
rium at the disk-corona boundary. 

The coefficients of viscosity and diffusivity are taken to 
be ay = 0.3 and o^d = 0.1 (as in R09). 

Boundary conditions. The boundary conditions at the in- 
ner boundary r = Rin are the following: The frozen-in con- 
dition is applied to the poloidal component Bp of the field, 
such that Br is fixed while Be and B^p obey "free" bound- 
ary conditions, dBe/dr — and dBd^/dr — 0. The density, 
pressure, and entropy also have free boundary conditions, 
d{...)/dr — 0. The velocity components are calculated using 
free boundary conditions. Then, the velocity vector is adjusted 
to be parallel to the magnetic field vector in the coordinate 
system rotating with a star. Matter always flows inward at the 
star's surface. Outflow of a wind from the stellar surface is not 
considered in this work. The boundary conditions at the exter- 
nal boundary r — Rout in the coronal region < < Odi and 
Od2 < < TV are free for all hydrodynamic variables. Here, Odi 
corresponds to the top surface of the and Od2 to the bottom 
surface. We prevent matter from flowing into the simulation 
region from this part of the boundary. We solve the transport 
equation for the flux function ^ so that the magnetic flux 
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Figure 8. Matter fluxes through a hemispherical surface of radius r = 10 above the equatorial plane (top set of curves) and below the equatorial 
plane (bottom set of curves) for cases with jld = 10 and different quadrupole moments of p.q = 20, 40, 60. 
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Table 1. Reference values for different types of stars. We choose the mass M, radius R^, equatorial dipole magnetic field Bd~), of the star and 
derive the other reference values for the case Jld = 10 (see §2.3). To apply the simulation results to a particular star one needs to multiply the 
dimensionless values from the plots by the reference values from this table. 



flows out of the region together with matter. If the matter has 
a tendency to flow back in, then we fix ^. In the disk region, 
restricted by two values of {Odi < < 9d2), we fix the den- 
sity at p = pd, and establish a slightly sub-Keplerian velocity, 
= HiQ(rd), where n = 1 — 0.003 so that matter flows into 
the simulation region through the boundary. The inflowing 
matter has a fixed magnetic flux which is very small because 

Rout ^ Rin- 

The system of MHD equations ( [1|4| ) was integrated nu- 
merically using the Godunov-type numerical scheme (see Ap- 
pendix of R09). The simulations were done in the region 
Rin < r < Rout, < < n. The grid is uniform in the 0- 
direction. The size steps in the radial direction were chosen so 
that the poloidal-plane cells were curvilinear rectangles with 
approximately equal sides. A typical region for investigation 
of asymmetric winds was 1 < r < 40, with grid resolution 
NrXNe — 104 X 80 cells. The simulation domain has 13 sHces 
in the radial direction and 10 sHces in 6 direction. Each sim- 
ulation run takes 4 — 12 days on 130 processors of the NASA 
high-performance facilities. Cases with stronger quadrupole 
component require longer simulations. The simulation time 
increases with increase of the quadrupole moment, and hence 



the longest simulations are those at jlq — 60. Test runs were 
also performed at the lower grid of 80 x 60 and higher grid of 
160 X 120 which show similar results (approximately the same 
matter flux onto the star and into the winds), though the lat- 
ter grids requires more computer resources. Simulations with 
the grid of 80 x 60 were used for a number of exploratory 
runs. 

3 RESULTS 

3.1 Properties of one-sided outflows 

Here, we chose one case with intermediate parameters: fid — 
10, pq — 20, call it "the main case", and show it in greater 
detail compared with other cases. Fig |2] shows formation of 
asymmetric outflows. One can see that initially, at t = 5 — 8, 
matter start to accrete onto the star along the shortest path, 
towards the quadrupolar belt, which is below the equator. 
Later, at t 15, the disk matter diffused through the external 
closed field lines of the top and bottom parts of the magneto- 
sphere, these lines inflated, and the conical-type wind start to 
blow along these field lines. Most of matter outflows from the 
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top side of the disk, where the magnetic field is stronger. Out- 
flows are episodic but quasi-steady on average. The episodic 
nature is connected with accumulation of matter near the 
closed magnetosphere, diffusion through closed field lines, 
inflation of these lines, outburst to the wind and enhanced 
accretion onto the star (see also Goodson & Winglee 1997; 
U06; R09). Figjsjshows larger view of typical, well-developed 
outflows at t = 50. One can see that outflows are more pow- 
erful on the top part, and much less powerful on the bottom. 

The mechanism of the outflow formation is similar to 
that in cases of conical winds (R09) : The magnetic flux of the 
star is pushed towards the star and compressed by the disk 
so that the poloidal field lines are always inclined relative 
to the disk. In addition, the rotating magnetosphere acts to 
make the field lines to rotate with a super-Keplerian velocity. 
These field lines thread the corona which rotates more slowly. 
Consequently there is a strong winding-up of the field lines 
just above the disk. Hence a strong magnetic force appears, 
Fm — — V(r^5|), which drives matter down the gradient of 
the magnetic pressure (Lovelace, Berk & Contopoulos 1991). 
Fig |4| shows the poloidal current through a disk of radius r, 
Ip oc rB(p. One can see that matter flows to the conical- type 
winds and is driven by the magnetic force. Ohmic heating is 
included in our code, but it does not have a significant role in 
the jet launching. Figure [s] shows the initial distribution of the 
specific entropy and that at t = 50. One can see that the low- 
entropy (cold) matter from the disk flows to the winds above 
and below the disk. On the top side of the disk this cold matter 
pushes the hot coronal gas towards the axis. On the bottom 
side of the disk, the hot coronal gas is pushed away and is 
gradually replaced by cold gas from the disk. 



3.3 Comparison of magnetic moments 

Below, we discuss the relative strengths of the dipole and 
quadrupole components of the field. For this it is useful to 
consider the dipole and quadrupole field components in the 



equatorial plane, Bdz 
radius at which \Bdz \ — \Bq 



-[idlr-^ and B^r = -3/ig/4r^. The 



3 

4 /id 



4 /id 2 /id 



(9) 



where we took into account that i?o = 2R^. We can esti- 
mate this radius for all cases above where /id = 10 and 
quadrupolar moments, \iq — \, 10, 20, 30, 40, 60. We obtain: 
= (0.15, 1.5, 3.0, 4.5, 6.0, 9.0)i?^. Hence, the quadrupole 
component is dynamically important in all cases, except /ig = 
1, and is expected to influence, e.g., the matter flow around 
the star at r > Veq (Long et al. 2009, 2010; Romanova et al. 
2009). However, we see that the quadrupole component de- 
termines the "one-sideness" and the direction of the matter 
outflow even for very small values, such as of /ig = 1. 



3.4 S)nTimetric Outflows for a pure quadrupole field 

We performed simulations of outflows in the case of a pure 
quadrupole field with /ig = 40, 60 (/id = 0). Figure [TT] shows 
simulation results at t = 50 for case with /ig = 60. In this case 



the outflows are almost symmetric about the equatorial plane. 
However, in a test case with very small dipole magnetic field, 
/i^ = 1, we found that the flow is one-sided and matter flows 
towards the direction of the larger intrinsic magnetic field. 



3.2 One-sided outflows for different quadrupole 
moments 

In this section we compare asymmetric outflows obtained in a 
set of simulations where we fixed the dipole component of the 
field at /id = 10 and varied the quadrupolar component from 
very small up to very large values: /xg = 1, 10, 20, 30, 40, 60. 
Figure |6] shows an initial magnetic field distribution in a num- 
ber of cases. Figure [6] shows that at /ig = 1 the quadrupole 
component is very small and the magnetic field is almost pure 
dipole field, while for /ig = 40, 60 the quadrupole field dom- 
inates. Figure [t] shows asymmetric outflows at t = 50. One 
can see that in all cases most of the matter outflows above the 
disk, to the side where the intrinsic magnetic field is stronger. 
At the same time there is a much weaker matter flux from the 
other side of the disk. 

We calculated the matter fluxes to outflows through a 
spherical surface of radius r = 10 

M = . , (8) 

where is the surface area element directed outward. Figure 
[s] shows matter fluxes above and below the disk for cases with 
different /ig. One can see that the main, upward outflows, 
have similar matter fluxes for /ig = 40 and /ig = 60, while at 
lower quadrupole moment /ig = 20 matter flux is higher. The 
matter fluxes of the downward outflows (bottom set of curves 
in Fig. |8| are similar for all three cases and are 4 — 8 times 
smaller than the main, upward outflows. 



3.5 "Flip-flop" of Outflows for a pure dipole field 

For a pure dipole field, the outflows behaved in an unexpected 
way. In past simulations the MHD equations were solved only 
in the upper half space and it was assumed that the full flow 
was given by reflecting the top flow about the equatorial 
plane (U06; R09). However, in the present work we calculate 
the disk and outflows in the entire space above and below 
the equatorial plane. We discover that during a brief initial 
time, t < 15, the outflows are symmetric about the equatorial 
plane (see fig. [9] top left panels) . However, later the outflows 
become strongly asymmetric with the direction of the main 
flow downward, which is an opposite to that of cases with 
quadrupolar component. Later, the main outflow changed its 
direction and matter started to flow upward. Later, the direc- 
tion changed again (see fig. [9| . The matter fluxes calculated 
through hemispherical surfaces with radius r = 10 above and 
below the disk also show this "flip-flop" behavior. The time- 
scale between the reversal of the direction of outflows is about 
30 rotation periods. For CTTSs this is about 30 days. 



4 CONCLUSIONS 

We performed axisymmetric MHD simulations of disk ac- 
cretion onto rotating magnetized stars. The star's intrinsic 
magnetic field was assumed to consist of a superposition of 
aligned dipole and quadrupole components. This field con- 
figuration is not in general symmetric about the equatorial 
plane. Thus the calculations must be done in both the upper 
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Figure 9. "Flip-flop"of outflows in case of a pure dipole field (/ig = 0, jld = 10. Color background shows the matter flux distribution, and lines 
are the magnetic field lines. 
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Figure 10. Matter fluxes calculated through the spherical radius r = 10 above the equatorial plane (solid lines) and below the equatorial plane 
(dashed lines) for the case of pure dipole field, jld = 10 and jlq = 0. 



and lower half-spaces. The star was assumed to be rapidly ro- 
tating with the magnetosphere in the propeller regime (U06). 
The ratio of the turbulent viscosity to the turbulent magnetic 
diffusivity in the disk was considered to be larger than unity 
so that the magnetic field threading the disk is advected in- 
ward (U06, R09). The main findings are the following: 

1. For cases with both dipole and quadrupole compo- 
nents, a one-sided conical wind forms and persistently blows 
in one direction. The favored direction is that of larger mag- 
nitude intrinsic axial magnetic field. Much weaker outflows 
form on the opposite side of the disk. 



2. For the case of a pure dipole field, the outflows are 
also one-sided, but the outflow direction alternates or "flip- 
flops" on a time-scale of about 30d for a T Tauri star. If the 
quadrupole component is small, jlq — I, the behavior is simi- 
lar to the case of /ig = 0. 

3. For the case of a pure quadrupole field, symmetric out- 
flows form. However, the presence of even very small dipole 
component leads to one-sided outflows. 

Note, that outflows from a T Tauri star may change direc- 
tion due to variations of the star's magnetic field on a time- 
scale of months (Smirnov et al. 2004). If the complex field 



© 0000 RAS, MNRAS 000, 000-000 





of a star is determined by dynamo processes inside the star, 
then one or another hemisphere may have the stronger axial 
magnetic field which determines the direction of the outflow. 

An intrinsic asymmetry of the outflows or jets will give 
a net force on the protostar by analogy with the previously 
analyzed case of asymmetric jets from magnetized black-hole 
disks (Wang et al. 1992; Tsygan 2007; Kornreich & Lovelace 
2008). If the outflow is one-sided for a long time T with 
asymptotic velocity Vj and mass flux Mj, then the velocity 
imparted to the star is A^;* = MjVjT/M^. For example, for 
a T Tauri star with Mj = lO~^M0yr~\ Vj = 3 x lO'^cm s~\ 
T 10^ yr, and M* Mq, we find A^;* = 3 x lO^cm s"^ 
which is probably undetectable. For the case of an almost pure 
dipole field the frequent "flip-flops" of the outflow direction 
will cause the star to random walk but the net displacement 
and velocity are very small. 

In addition to the intrinsic stellar magnetic field con- 
sidered here, the accretion disk can advect inward exter- 
nal (e.g., interstellar) magnetic flux because of the disk's 
highly conducting (non-turbulent) surface layers (Bisnovatyi- 
Kogan & Lovelace 2007; Rothstein & Lovelace 2008; Lovelace, 
Rothstein, & Bisnovatyi-Kogan 2009). The combination of 
the advected field and the intrinsic field of the star can 
give rise to a complex field structure near the star's magne- 
topause which produces asymmetric or one sided outflows. 
The field of the star may be dynamo generated with a complex 
time-dependent structure (e.g., von Rekowski & Brandenburg 
2006) . In the case of disk accretion to a black hole, a large- 
scale asymmetric magnetic field close to the black hole can 
arise from advection of external flux due to the conducting 
surface layers of the disk or it may arise from dynamo pro- 
cesses in the disk which generate both dipole and quadrupole 
field components (e.g., Pariev & Colgate 2007; Pariev, Col- 
gate, & Finn 2007). Dynamo processes may also be important 
in the disks of accreting stars. 
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